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The baryon-strangeness correlation in the hadronization of the quark matter is studied within the quark com¬ 
bination mechanism. We calculate the correlation coefficient Cbs = “3((BS) - {B){S))/{(S^) — {S)^) of initial 
hadrons produced from the deconfined free quark system with Cgj = 1. The competition of the production of 
baryons against that of mesons is the key dynamics that is most relevant to the change of baryon-strangeness 
correlation during system hadronization. Results of quark combination under the Poisson statistics agree with 
the statistical model predictions for a hadron resonance gas at vanishing chemical potential but dilfer from those 
at relatively large chemical potentials. Results beyond Poisson statistics are also obtained and are compared with 
calculations of Lattice QCD in the phase boundary, giving the best agreement at temperature T = 163 MeV. We 
predict the dependence of the Cbs of hadron system on the baryon chemical potential and strangeness. These 
predictions are expected to be tested by the future Lattice QCD calculations at non-zero chemical potentials 
and/or by the Beam Energy Scan experiment of STAR Collaboration at RHIC. 

PACS numbers: 25.75.Gz, 25.75.Nq 


I. INTRODUCTION 


Baryon-strangeness correlation is an effective diagnostic 
tool for the relevant degrees of freedom of the hot nuclear mat¬ 
ter produced in relativistic heavy ion collisions . It has been 
extensively investigated by various phenomenological models 
of high energy collisions IJ-Ql and first-principle calculations 
in Lattice QCD The correlation is usually quantified by 
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Here, angle brackets denote the event/ensemble average. In 
the second step, strangeness neutrality (S) = 0 is applied 
to relativistic heavy ion collisions. For a deconfined sys¬ 
tem consisting of free quarks and antiquarks, - 1 be¬ 
cause strangeness is carried only by the strange (anti-)quarks 
which carry the baryon number in strict proportion to their 
strangeness with the coefficient -1/3, i.e., Bs = In 

contrast, the relation between baryon number and strangeness 
in a hadron system is multiple, e.g., baryon number 1 for a 
strange baryon but 0 for a strange meson. Cbs of hadron sys¬ 
tem is usually smaller than one (at low baryon number den¬ 
sity) due to the fact that most strange quarks, at hadroniza¬ 
tion, will come into mesons instead of baryons which unlock 
the intimate baryon-strangeness correlation existed previously 
in quarks. Statistical model estimation of Cbs for a hadron 
resonance gas at zero baryon number density is about 0.66 
111] . Calculations of Lattice QCD also show that the Cbs of 
the strong-interacting system at temperature above the phase 
transition temperature Tc tends to one while near and below 
Tc the system Cbs decreases rapidly @,|^. 

Hadronization refers to the process of the formation 
of hadrons out of quarks and/or gluons, accompanying 
with the change of the correlation between baryon num¬ 
ber and strangeness. Phenomenological models describing 
the hadronization should reproduce this change of baryon- 
strangeness correlation due to the transformation of the de¬ 
grees of freedom in the system. Quark combination is one 


of the effective mechanisms for the hadronization of the hot 
quark matter produced in relativistic heavy ion collisions and 
has explained lots of experimental phenomena in heavy ion 
collisions at SPS, RHIC and recent LHC, see, e.g. RefsH^ 
[lsIi . The change of baryon-strangeness correlation during 
quark combination hadronization is intuitive. A strange 
quark(antiquark) combines with a light antiquark(quark) to 
form a strange meson, which completely destroys the origi¬ 
nal baryon-strange correlation carried by strange (anti)quark, 
i.e., B^sq) =0x5 (sq), but inherits the strangeness lossless. 
Occasionally, a strange quark combines with a strange anti¬ 
quark to form a hidden-strange meson, which inherits nothing 
from original quarks. On the contrary, in the baryon forma¬ 
tion a strange quark combines with two light quarks to form 
a baryon, which alters the baryon-strangeness correlation co¬ 
efficient B^sqq) = -S{sqq)- The Combination of two strange 
quarks with one light quark also alters the correlation coef- 
hcient B^^sq) - -jS^ssq)- Obviously, the Cbs of the system 
depends on the relative proportion of the produced baryons to 
mesons. 

In this paper, we study the change of system Cbs caused 
by hadronization in quark combination mechanism (QCM). 
We discuss in details how the dynamics of baryon-meson pro¬ 
duction competition at hadronization to dominate the Cbs of 
hadron system. In addition, we study the dependence of Cbs 
on the strangeness content and the baryon number density of 
the system, which is related to the hot quark system produced 
in relativistic heavy ion collisions at different collisional ener¬ 
gies. We compare our results with the calculations of Lattice 
QCD 0,0] and the prediction of statistical model for a hadron 
resonance gas Q]. 

The paper is organized as follows. In Sec. II, we present a 
working model in QCM for the yields of hadrons. In Sec. Ill, 
using the model we explain the experimental data of the mid¬ 
rapidity yields of strange hadrons that are most relevant to 
Cbs calculations in central Pb-i-Pb collisions at s/Imf = 2.76 
TeV. In Sec. IV, we calculate the Cbs of the hadron system 
in which two assumptions are made for initially produced 
hadrons. One assumption is that different kinds of hadrons 
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are uncorrelated, which can be expected in the case that the 
quark system existed previously is made up of free quarks and 
antiquarks = !)■ In this paper, we only consider this kind 
of quark system since the effects of hadronization on system 
Cbs are already addressed properly. The other assumption is 
Poisson statistics for yield fluctuations of hadrons. In Sec. V, 
we present a different approach of calculating Cbs which is 
beyond Poisson statistics and can also reflect intuitively the 
essence of baryon-strangeness correlation. Summary and dis¬ 
cussion are finally given in Sec. VI. 


II. A WORKING MODEL OF HADRON YIELDS 


In this section, we present a working model of the yields of 
hadrons that are produced from the deconfined quark phase 
after hadronization in the framework of quark combination 
mechanism, for the facility of the subsequent correlation stud¬ 
ies. As discussed in introduction, baryon-strangeness correla- 
tion is strongly dependent on the relative production of 
baryons to mesons, which requires our model to well address 
this point. Therefore, we introduce the working model ac¬ 
cording to the following strategy: (1) give the global proper¬ 
ties for the global production of all mesons, all baryons and all 
antibaryons, in which the production competition of baryons 
against mesons is properly addressed. (2) on the basis of (1), 
give the yield formulas of various identified hadrons. 

In previous work ifl^ . we have studied the properties of 
the global production of all mesons, baryons and antibaryons, 
and obtained their yield formulas utilizing only the basic ideas 
of QCM. We start from a deconfined system consisting of Nq 
constituent quarks and Nq antiquarks. Here, gluons contribu¬ 
tion to the system on the threshold of hadronization is replaced 
by pairs of quark and antiquark. After hadronization, the sys¬ 
tem changes the basic degrees of the freedom to become the 
hadronic system and produces in average BiNq, Nq) baryons, 
B(Nq,Nq) antibaryons and M(Nq,Nq) mesons. Four proper¬ 
ties of hadron production from general principles are used to 
constrain the behavior of their yield formulas. 

(1) charge conjugation symmetry of the hadron yields, 

B(Nq,Nq)^B(Nq,Nq), 

M(Nq,Nq)^M{Nq,Nq). 

(2) unitary of the hadronization, i.e. production of mesons, 
baryons and antibaryons should exhaust all quarks and anti¬ 
quarks of the system existed previously 


M(Nq,Nq) + 3B(Nq,Nq)^Nq, 
M(Nq,Nq) + 3B(Nq,Nq)^Nq. 

(3) boundary condition, i.e., 

(b^O, B^ ^,M = 0 ifA^ = 0 
= M = 0 ifA^ = 0 


(4) linear response of meson and baryon yields to quark anti¬ 


quark numbers 


M(ANq,ANq) = AM(Nq,Nq), 
B(ANq,ANq)^AB(Nq,Nq), (5) 

B(ANq,ANq) = AB(Nq,Nq). 

Using these properties, we obtained in Ref. lfl^ the yield for¬ 
mulas of baryons, antibaryons and mesons. 
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Here, we have rewritten x — Nq + Nq which characterizes the 
bulk property of the system related to the system size or en¬ 
ergy and rewritten z = (Nq-Nq)lx which depicts the asymme¬ 
try between quarks and antiquarks in the system (|z| < 1) and 
is a measurement of the baryon number density of the system. 

The production competition of baryons against mesons is 
often quantified by the yield ratios RB/miz) = B(x, z)/M(x, z) 
and Rbim{Y) = B{x,z)IM{x,z), which are the function of only 
z. The factor a in Eq. (|6]l represents the degree of the baryon- 
meson competition by the relation a = ^ ■ Studies in 

Ref. im show the Rb/m{0) of value about 1/12, i.e., a 5, 
can well describe yield ratios of various baryons to mesons in 
heavy ion collisions at LHC energy. We note that such a com¬ 
petition as well as yield formulas Eq. (|6]l can be properly ad¬ 
dressed by a phenomenological combination rule in the quark 
combination model developed by Shandong group (SDQCM) 
ifT^ . In addition, using vested Rb/m( 0), we have successfully 
explained the yield ratios of various anti-hadrons to hadrons at 
non-zero z region in relativistic heavy ion collisions 
i.e., the data of these yield ratios at different collision energies 
and at different rapidities. We emphasis that through the de- 
pendence of Cg^ on the factor a or Rb/m{0) we can address the 
effects of hadronization on the baryon-strangeness correlation 
of the system. 

To study C^gg , we have to obtain the yield formulas of iden¬ 
tified hadron. Given the total yield of baryon, that of an- 
tibaryon and that of meson, the inclusive/averaged yields of 
identified hadrons Mi(qiq 2 ), Bjiqiq 2 qY) are calculated by their 
individual production weights. 


Nm, = PMiMix,z) ^ CM,Pqtq2.MMix,z), (7) 
Nbj = PbjB(x,z) ^ CBjPqiq2q,,BB(x,z). ( 8 ) 

Here, Pq,q 2 ,M denotes the probability that, as a meson is known 
to be produced, the flavor content of this meson is qiq 2 - Cm; 
further denotes the branch ratio of this meson with given flavor 
composition qiq 2 to a specific meson state M,-. CMiPqiq 2 ,M 
thus gives the emerging probability of a specific meson, Rm, , 
when a meson is known to be formed. Similarly, Pqiq 2 q 2 ,B 
denotes the probability that, as a baryon is known to produced, 
the flavor content of this baryon is q\q 2 q 3 , and Cb, denotes 
the branch ratio of this baryon with given flavor composition 
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qiqiqi to a specific baryon state Bj, and thus gives 

the emerging probability of a specihc baryon, Pg., when a 
baryon is known to be formed. A similar formula holds for 
antibaryons. 

The probability Pq^q 2 ,M is posteriorly evaluated by the pro¬ 
portion of qiq2 pairs in all quark-antiquark pairs in the sys¬ 
tem, Pq^q^.M - Nq^q^INqq, by Considering the fact that ev¬ 
ery quark/antiquark has both the probability of entering into 
a meson and the probability of entering a baryon/antibaryon. 
Here, Nq = Y^iNq. and Nq - YiNq^ are total number of 
quarks and antiquarks in the system, respectively. We have 
^q\qi - ^qx^qi Nqq - NqNq. Similarly, the probabil¬ 
ity Pq^q^q^^B IS evaluated by the proportion of q\qzq-i combina¬ 
tions to all three-quark combinations in the system, Pq^q 2 q,,B - 
NiterNq.q^qJNqqq. Nq^q,q^ is the number of qiqzqz, combina¬ 
tion which satisfies Nq^q^q^ - Nq^Nq^Nq^ for qx + qz + qz, 
= ^qA^q^ - for q\ = qi * q3 and Nq^q,q, = 

Nq^Nq^ - l)(Nq, - 2) for qi ^ qz = qz. Nqqq = - 

^)iNq - 2) is the number of all qqq combination in the sys¬ 
tem. In large quark number limit we can apply Nqqq » 
and Nq^q^q^ ^ Nq^Nq^Nqj always. Niter stands for the number 
of possible iterations of q\qzqz, which is 1, 3, and 6 for three 
identical flavor, two different flavor, and three different flavor 
cases, respectively. 

Cmj and/or Cb^ denote the branch ratio of a given flavor 
composition (e.g., us) to a specific hadron state (e.g., K^) un¬ 
der the condition that they are known form a hadron. In the 
case that when the ground state - 0“ and fo mesons and 
\ and I baryons are considered only, we have, for 
mesons, 

1/(1 -H Rvjp) for - 0“ mesons, 

Rvjpli^ + Rvjp) for - fo mesons, 

where Rvip represents the ratio of the - 1“ vector mesons 
to the - 0“ pseudoscalar mesons of the same flavor com¬ 
position; and for baryons, 

RoidH 1 + Roid) for 7^ = (1 /2)+ baryons, 

1/(1 -H Roid) for 7^ = (3/2)+ baryons, 

except that C\ - C^o = IRqidK^ +'^Roid), = 

1/(1 -H 2Ro/d), Ca++ = Ca- = Cq- = 1 . Here, Rq/d stands for 
the ratio of the 7^ = (1/2)+ octet to the 7^ = (3/2)+ decuplet 
baryons of the same flavor composition. The two parameters 
Rv/p and Rq/d can be determined using the data from different 
high energy reactions ifT^flsIl . and they are taken to be 3 and 
2 in this paper, respectively. 


III. YIELDS OF STRANGE HADRONS AT LHC 

With above working model, we can conveniently predict 
yields of various identified hadrons. As an illustration, we 
now explain the experimental data of mid-rapidity yields of 
strange hadrons K, A, H and in Pb-i-Pb collisions at 
s/snn - 2.76 TeV. These hadrons are most relevant to Cbs 
calculation. 


As applying yield formulas in Sec. II to the hnite rapid¬ 
ity window in heavy ion collisions, we note that Eqs. (|2]l 
and (|4]l5]l still hold generally and Eq. Q also approximately 
holds with good numerical accuracy due to the locality of the 
hadronization, and therefore no conceptual issues exist for our 
mid-rapidity predictions. In addition, our formulas in Sec. II 
give the averaged yields of hadrons at the given quark system 
with the fixed quark numbers while the experimental data of 
hadronic yields are event-averaged quantities. The produced 
quark system in heavy ion collisions at a given collision en¬ 
ergy is varied in size event-by-event and thus the number of 
quarks and that of antiquarks should follow a certain distri¬ 
bution around the averaged quark number {Nq.) and antiquark 
numbers (Nq.) (where i = u,d,s considered in this paper). 
Our final predictions of hadron yields should be the average 
over this distribution. In general such averages depend on the 
precise form of the distribution. Here, we approximate these 
averages by taking the corresponding values of the quantities 
at the event averages {Nq.) and {Nq^), i.e.. 


{NMt) = CM.PqtqtMM({x),{z)), (9) 

{Nb) = CB;Pqtqtq,,BB({x),{z)), ( 10 ) 


where {x) - {Nq)+{Nq), {x){z) - {Nq)-{Nq) and P is also cal¬ 
culated with the event averaged quark numbers. Such approx¬ 
imation can be expected in considering that properties Eqs. (|2] 
- nil in our derivation of hadron yields are also apparently sat¬ 
isfied even for the event-averaged quantities. 

Eor convenience, we use factors A., - {Ns)/{Nu) = 

{Ns)/{Nj) and associated A{ = <A^)/<A„) = {N,)/{Nd) to de¬ 
note the suppression of strange anti-quarks relative to light 
antiquarks and that of strange quarks to light quarks, respec¬ 
tively, and simply the yield formulas of hadrons. Here, isospin 
symmetry between (anti-)up quark and (anti-)down quark is 
applied. 

Considering the decay contribution from the short-lived 
resonances, we obtain the yields of final state hadrons 


= {Nh.) + ^ Br{hk ^ hj) {Nitt), ( 11 ) 

k 


where we use the superscript (/) to denote the results for the 
final hadrons to differentiate them from those for the directly 
produced hadrons. The data of the decay branch ratios are 
taken from the PDG ifl^ . 

With = 7' at LHC ({z) - 0 approximation), we obtain 
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yields of these strange hadrons in final state 
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Here, we have only taken into account the strong and electric- ^ g-3 

magnetic (S&EM) decays of short-lived resonances. 

Absolute yields of these hadrons are dependent on the total 
quark number (x) of the system by M((x), 0) = 2{x)/5 and 
B({x},0) - (x)/30 with fixed/? b/m( 0) = 1/12. To eliminate 
this (x) dependence, we consider the relative production of 
these hadrons to pions, i.e., the yield ratios of these hadrons 
to pions, which finally rely only on the strangeness of the sys- q-4 

tern, beside of the baryon-meson competition in hadroniza- 
tion. Since the decay contributions to pion are complex, we 
here give directly the numerical result of the yield of i.e.. 



0 50 100 150 200 250 300 350 400 



- 0.213{x) under S&EM decays, instead of the detailed 
compositions like Eqs. (I12II15I I. In Fig.[T] we use Eqs. (I12II15I I 
to explain the experimental data of the mid-rapidity yield ra¬ 
tios In^, AItt^, E~In^ and Dr/n^ in central Pb-i-Pb colli¬ 
sions at ^jsNN - 2.76 TeV ||351 - It^ . To incorporate the change 
of Tj at different collision centralities, a varied strangeness 
AsiNpart) - (0.43 + 0.02)/(1 -h lO.SA'J,:^) is used. We can see 
that the hierarchy properties in yields of these strange hadrons 
and their Apart dependence can be systematically described 
by our formulas Eqs. (I12II151 I. Here, we would like to empha¬ 
sis that the yield difference between K and those hyperons is 
mainly due to the baryon-meson competition at hadronization 
while the hierarchy structures among A, and D are mainly 
strangeness relevant. 


FIG. f: (Color online) Yield ratios Ain*, "Sr In* and In* 

at mid-rapidity as the function of nuclear participants Npan in central 
Pb+Pb collisions at yIsNN - 2.76 TeV. Symbols are experimental 
data from Refs. 13514^ . Shadow regions are our results with a Ap„„- 
dependent strangeness AANpan) = (0.43 ± 0.02)/(l -l- lO.SAp/;/). 


IV. Cbs of hadrons under POISSON FLUCTUATIONS 

For initial hadrons produced by the hadronization of de- 
confined quark system, the baryon number of the system is 
B - Yja Qa.sAa and the strangeness S - Qa,sAfa, where 
the species a has baryon number Qa,B and strangeness Qa,s ■ 
By definition Eq. ([T]i, the Cbs of hadrons is 


Applying our yield formulas to other hadrons at LHC and 
those at RHIC energies with nonzero baryon number densi¬ 
ties, we also find a good agreement with available experimen¬ 
tal data. This is not surprise. In fact, QCM has already shown 
its effectiveness in explaining the data of hadronic yields and 
longitudinal rapidity distributions in relativistic heavy ion col¬ 
lisions at different collisional energies ifT^flsl 12^241] . The re¬ 
lated low-pj- issues of QCM such as entropy conservation and 
pion production have been properly addressed ifl^ 1^ l25l - 
1^ . There are also many successful applications of QCM on 
correlation studies, e.g., multi-hadron yield correlations iS- 
[HHl, baryon-meson correlated emission ll^l^ as well as 
the charge balance function 


Ya,/3 Qa,S Qp,S Ca/j 

where the covariance Ca/s - {NaN^} - {Na}{N/j} describes 
the correlation between hadron a and hadron yS. A general 
calculation of in QCM is still unavailable in the current 
progress of hadronization phenomenology since there still has 
lots of unsolved dynamics in hadronization due to its non- 
perturbative feature. On all “on market” combination mod¬ 
els, there are a few ones that can give the calculation of Coy? 
with their own specific model details/assumptions. Such spe¬ 
cific calculations are not the purpose of this paper since we 
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intend to analyze the hadronization effects in a general and 
transparent way. Here, we consider a simple case that after 
hadronization different kinds of produced hadrons are uncor¬ 
related, i.e. Cajs = Sajjcr^- This can be expected if the quark 
system existed previously is made up of free quarks and an¬ 
tiquarks (i.e., = !)■ We note that, above Tc, the Crs of 

strong-interacting system calculated by Lattice QCD in¬ 
deed tends to one and the off-diagonal flavor susceptibilities 
also tend to be relatively small. We only consider this kind of 
quark system in this paper because the effects of hadronization 
on system Crs are already addressed clearly . 

In this section, we first calculate the Cbs of the initial 
hadron system, assuming the Poisson statistics » (AL) for 
the yield distribution of identified hadrons. Now we have 

_ ^ YjaQa,BQa,s{Na) 

With Eqs. (|9]l and ( fTOl i. we are ready to calculate Crs of 
hadrons via Eq. (El), 

{BS) = -{{Nk) + {N^.,) + {N^..,)]-2{{N^..-) + {N^,..-)] 


about 5% as As increases from 0.3 (the rough value in pp re¬ 
actions) to 0.7 (almost maximum value occurred in heavy ion 
collisions). In addition, we see that of initial hadron sys¬ 
tem is obviously small than one (the value of ideal quark sys¬ 
tem) because the produced strange mesons significantly out¬ 
number the strange baryons in the current baryon-meson com¬ 
petition Rb/m{0) - 1/12. Eig. 2 (b) shows the dependence 
of on baryon-meson competition factor /?b/m( 0), as the 
strangeness is taken to the saturated value 0.43 in relativis¬ 
tic heavy ion collisions. Here, the saturated dj is extracted 
from the fit of yields of strange hadrons in Sec. Ill and in our 
previous work IT^ . and we note that this value is consistent 
with calculations of Lattice QCD at Tc iSH^. Clearly, we see 
that the increase of Rbim{Q) will enhance the yields of baryons 
against mesons and the increases rapidly. With preferred 
values RB/uiO) - 1/12 and d^ = 0.43, is about 0.65, 
much smaller than = 1. To reach the unit correlation for 
hadrons, an extremely high baryon-meson competition factor 
Rbim(0) ~ 1/5 is needed, which is completely unable to re¬ 
produce yields of strange mesons and baryons. 



( 20 ) 

ih') 

which gives the dependence of on baryon-meson compe¬ 
tition factor Rbim(0), strangeness dj and the baryon number 
density of the system. 

We firstly consider the situation of zero baryon number den- 
sity {z} = 0 to study the dependence of on Rbim(0) and 
As. With d/ = As and Rb/m(0) = Rb/m(0)’ we get 

rW + As)Rb/m{0) 

(3As + 2)Rb/m{0) + 2I3' ^ 

Fig. 2 (a) shows the dependence of C® on the strangeness dj, 
as the Rb/m(Q) is taken to 1/12. One can see that is insen- 
sitive to the change of the strangeness, i.e. only increases 


FIG. 2: (Color online) (a) The dependence of Cgj on strangeness 
suppression factor d,; (b) on baryon-meson competition factor 
Rb/m(0). 

Subsequently, we study the dependence of on the 
baryon number density of the system. In previous discussions, 
we use the quark-antiquark asymmetry {z} to characterize the 
baryon number density of the system. To compare our results 
with existed predictions of statistical models, we alternatively 
use the chemical potential pB which relates (z) via 

2 sinh(2^) 

<^> = - u.Mbs ; r MSS ’ (22) 

2cosh(^) -H As exp(-^) 
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under the assumption of Boltzmann distribution for thermal- 
ized quarks and antiquarks. Here, strangeness neutrality Ns - 
Ns is applied. Tc is the temperature of the quark system at 
hadronization. 

Fig. [3 shows our predictions of the of initial hadrons 
as the function of Here, we have taken into account the 
fact that both Tc and are varied with the /jg for the hot 
quark matter produced in relativistic heavy ion collision. For 
the fig dependence of Tc, we apply the calculation of Lat¬ 
tice QCD by G. Endrodi et al i.e., Tc(p.g) = F£.o(l - 
0.0089^(g/r^g) with Tco being 0.169 GeV, the transition tem¬ 
perature at fig = 0. For the /ug dependence of As, we use our 
previous extractions at mid-rapidity in heavy ion collisions at 
different collision energies ifo HM . i.e., dj = (0.43, 0.43, 
0.44, 0.44, 0.48, 0.5, 0.57, 0.8, 0.7) as = (2760, 200, 

130, 62.4, 17.3, 12.3, 8.7, 7.6, 6.3) GeV, and convert the en¬ 
ergy dependence into fig dependence by the formula - 

(1.308GeV//rB - l)/0.273 GeV in Refs.® IHl]. The solid 
line with filled square symbols are our result. We also calcu- 
late Cgg at hxed d^ = 0.43 (dotted line with down-triangles), 
at fixed Tc - 0.169 GeV (dashed line with up-triangles), and 
at both fixed Tc - 0.169 GeV and fixed d^ = 0.43 (dashed line 
with open circles), respectively. 



FIG. 3: (Color online) The dependence of Cbs on the chemical 
potential fig of the system. The lines with filled squares, down- 
triangles, up-triangles and open circles are our results at both varied 
dj and Tc, at constant dj = 0.43, at constant Tc = 0.169 GeV, and 
at both constant d., = 0.43 and Tc = 0.169 GeV, respectively. The 
solid line with open cross symbols is our results using the chemical 
freeze-out temperature. They are compared with the prediction of 
statistical model for hadron resonance gas (T|], the dashed line with 
stars. 

A striking behavior of of initial hadrons is that it in¬ 
creases with the increasing fig and in the large fig region 
ifiB k, 0.3 GeV) it surpasses the unity and become higher at 
larger fig. This is in sharp contrary to the free quark system 
existed previously where the cotTelation coefficient remains 
strictly unity at all temperatures and chemical potentials. This 
is because that, as the fig increases, the relative production 
of baryons to mesons, i.e., Rg/miz) in our approach, becomes 


large, and then weight of this item in Cgg formula Eq. (l20l i 
increases and thus Cgg increases correspondingly. Compar¬ 
ing the result of Cgg at fixed Tc with those with varied Tc, we 
see that the change(decrease) of Tc at large chemical poten¬ 
tial increases the Cgg of the system about several percentages. 
Comparing the result of Cgg at fixed dj with those with var¬ 
ied dj, we find that the change of the strangeness also slightly 
influences the Cgg of the system. 

We also comp^e our results with the early prediction of 
V. Koch et al (Tl] (the dotted line with star symbols) for 
a hadron resonance gas in the statistical treatment. A fig 
dependent chemical freeze-out temperature, TchemicaiiRB) ~ 
0.17 - 0.13/i| - 0.05yu^ in Ref. id^l . is used in their predic¬ 
tion, which falls obviously below the transition temperature 
we used above at intermediate and large fig. For comparison, 
we also calculate the Cgg of initial hadrons using their chemi¬ 
cal freeze-out temperature and the result is presented in Fig. [2 
as the solid line with open cross symbols. We find that our 
result is close to the prediction of statistical model at small 
fig but at intermediate and large fig region our result is about 
15-20% larger than that of statistical model. We note that re¬ 
cently F. Becattini et al ll43l] have reconstructed the original 
hadro-chemical equilibrium temperature after considering the 
effects of final hadron expansion phase and the results have 
closely followed the phase transition boundary line predicted 
by Lattice QCD j^ . 

V Cbs of HADRONS BEYOND POISSON STATISTICS 

In above calculations of , Poisson statistics for hadronic 
yields is applied as an open approximation. Compared the ob¬ 
tained at zero fig with calculations of Lattice QCD in the 
phase boundary |j^, we find our result conform to Lattice cal¬ 
culations at r = 150 + 3 MeV. This value is lower than the 
usually expected deconfined temperature (around 165 MeV) 
and is also lower than the chemical freeze-out tempera¬ 
ture extracted from hadronic yields in statistical models which 
is about 166 MeV l^ldd]] . Such an underestimation may be 
relevant to the applied Poisson statistics approximation since 
this application implies that a further sufficient evolution of 
hadron phase dominated by inelastic scatterings is probably 
needed. Therefore, in this section, we calculate the of 
initial hadron system by another approach which is indepen¬ 
dent of the Poisson statistics assumption and also reflects in¬ 
tuitively the essence of the baryon-strangeness cotTelation. 

To probe the intrinsic cotTelation between the baryon num¬ 
ber and strangeness for a hadron system, we suppose that the 
system has a change by stochastically emitting a small amount 
of strange hadrons and then the system strangeness changes 
and the baryon number also changes Regarding 

the baryon number of the system as the function of the 
strangeness S‘'^\ the change of baryon number due to the per¬ 
turbation of strangeness can be calculated by 

d(^) 

= (^) 55''’) + Q((55<'’))2), (23) 

(fl) 

where Bf is the baryon number carried by strange hadrons. 
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To second order fluctuation of strangeness we get 


bined into hadrons after hadronization. Then, we have 




(24) 


The partial derivative is evaluated at the event average values 
of hadron yields {Nh)- Note that here we don’t consider the 
correlations between light and strange hadrons induced by the 
interactions among hadrons, i.e., we consider an ideal hadron 
gas system after hadronization. 

The same philosophy can be applied to the quark system 
before hadronization. Considering that the quark system is 
made up of free quarks and antiquarks with three flavors, 
the strangeness of the system = —{Ns - Ns) and the 
baryon number carried by these strange quarks and antiquarks 
- \{Ns - Ns) and we get = 1 exactly. The difference 
of Cbs between quarks and initial hadrons can be illustrated 
by the different slopes of two phases in B-S plane, as shown in 
schematic Fig. |4] The cross point between two phases stands 
for the hadronization that changes the basic degrees of free¬ 
dom of the system and also stands for the global charge con¬ 
servations during hadronization. 



FIG. 4: (Color online) Schematic picture for the transition of ideal 
quark phase to hadron phase by hadronization in S - S plane. 


With yield formulas Eqs. (|7]i and ([8]l, we have 


fiW = Na H- N^±fi + -h Nho,- -h Nh.o,- -h Na- 

— anti — hyperons (25) 

(>Nuds + 'iNuus + ^Ndds + ^Nuss + 'iNdss + Nsss 


Nqqq 

^Nms + 3Nuus + ^N^js + ^Niiss + ^N^ss + Nsss 

Nqqq 


-B{x,z) 


B{x, z). 


and the strangeness of initial hadrons 5^^' = YjaQa,sNa - 
-{Ns-Ns) because all strange quarks and antiquarks are com- 


Cfs - 3 X 


24B(1,<z» 


24.B(1,<z» 


(2 + T')3(1+<z» (2 + T,)3(1-<z» 

12T -H 6T? -H 64'? 


(2 + 4)3 




a-l 


[(1 - - 2a(z) -!) + (!+ (z))"^*] 

3[(l+<z»“-(l-<z»''f 

12/lv -t" 64 + 6/1^ . 

1(1 _ 


(2 + 4,)3 

[d + + 2a{z} - 1) + (1 - 


3[(1 + (z))" - (1 - <z»T 


(26) 


for initial hadron system. This result shows a complex de¬ 
pendence on baryon-meson competition factor /?b/m( 0) by the 
factor a = sRgjsiiO) ^'^'angeness 4^ and the baryon number 
density of the system by z. 

In the case of zero baryon number density z = 0, we have 
4' = 4, and B(l, 0) = B(l, 0) = l/6a and simplify Eq. (l26l l 


Mh) ^ , 6Rb/m{0) - 1 

(2 + 44(1+3/?b/m(0))' 


(27) 


Fig.E] (a) shows the dependence of on the strangeness 
4j, as the Rb/m{0) is taken to 1/12. We see that of 
initial hadrons in the current approach is also insensitive to 
the strangeness of the system, but the result is about 20% 
larger than the previous prediction in Eig. |2] (a) in Poisson 
fluctuation. Eig. |5](b) shows the dependence of Cbs on the 
baryon-meson competition factor/? b/m( 0). We see thatC^'^ is 
strongly dependent on the Rb/m{0), which is also qualitatively 
consistent with the result in Eig.|2](b). 

The qualitative consistency between Eig.|5]and Eig.|2]can be 
naturally understood via the definition Cbs - -^{BS}/{S^). 
Because both the nominator {BS} and denominator (S^) are 
dependent on the strangeness mainly by the single-strange 
hadrons, the strangeness dependence are thus partially offset 
in their ratio and we observe a weak dependence of on 
4j. On the other hand, since the nominator (BS) receives the 
larger contributions from multi-strange baryons than the de¬ 
nominator (5 3) which is always dominated by single-strange 
mesons, we would observe a slightly increase of with the 
increasing 4^. In addition, the increase of Rbim{0) enhances 
the production weights of (strange) baryons against mesons 
and contributes largely to the nominator {BS ) and we observe 
a rapid increase of with the increasing Rb/m{0). 

Since results in current approach hold more generally than 
those in Sec. IV under Poisson statistics, we can compare our 
result ^olid line in Pig. |5](b)) with the calculation of Lattice 
QCD 0 (box area) in the phase transition boundary. With 
the critical temperature Tc - 0.16 - 0.17 GeV, Lattice QCD 
predicts the Cbs of the system is about 0.74 - 0.85, which 
has a striking overlap with our results as Rbim{0) is within 
[1/9.5,1/13.5]. In the preferred situation of Rb/m{0) =1/12 
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FIG. 6: (Color online) The dependence of Cbs on the chemical 
potential jig of the system. The filled squares, up-triangles, open 
squares and open circles are our results at both varied T., and Tc, at 
constant T., = 0.43, at constant Tc = 0.166 GeV, and at both constant 
T, = 0.43 and Tc = 0.166 GeV, respectively. They are compared with 
the prediction of statistic method for a hadron resonance gas (I|]. 


FIG. 5: (Color online) The dependence of Cbs on strangeness sup¬ 
pression factor /Ij, (a), and on baryon-meson competition factor 
Rb/m(0), (b). Predictions of Lattice QCD are from Ref.j^. 


and Tj = 0.43, our calculated Cbs is 0.78, which best con¬ 
forms to the Lattice QCD calculations Cbs — 0.77 - 0.80 at 
Tc = 0.163 GeV within its theoretical uncertainty fl. Such a 
agreement with Lattice QCD calculations suggests that quark 
combination mechanism reflect some basic dynamics of the 
hadronization process of the strong-interacting quark system 
and thus it is indispensable for the hadronization phenomenol¬ 
ogy- 

We also predict the dependence of Cbs on the baryon num¬ 
ber density of the system using the chemical potential The 
results are shown in Fig. |6] Here, the transformation of vari¬ 
able z to fj.B is by Eq. (l22ll and the fiB dependence of Tc and Tj 
are same as those in Sec. IV. The results of statistical model 
from Ref. Q] at different /jb are also presented as the com¬ 
parison. The solid line with filled squares are our results at 
varied /jb and Tc. It is shown that they are about 15% higher 
than the predictions of statistical model, i.e., dashed line with 
stars, in the full region. Comparing our results at varied 
fiB and Tc with those at fixed fiB and/or Tc, we find that de¬ 
crease of Tc at large chemical potential will slightly in¬ 
crease the Cbs of the system and that the change of system 
strangeness at different ^b slightly influences the Cbs of the 
system. These effects of and on the Cbs of the system 
are similar to those calculated in Sec. IV. In addition, we see 
that results of small region calculated here are higher than 
those in Fig.[3]in Sec. IV but the results at large are close 
to the previous ones. We also note that another difference be¬ 
tween the current calculations of Cbs and the previous ones 
in Fig. |3]is the curvature of Cbs as the function of The 
Cbs in Fig.|3]is the concave function at intermediate and large 


jiB while the current Cbs in Fig.|6]is almost linear function at 
intermediate and large By noticing that the predictions of 
statistical model in this region is a slightly convex function, 
these different predictions may be identified in future by cal¬ 
culating the higher-order susceptibilities for the conservative 
charges of the system. The other difference between QCM 
prediction and statistical model is the position that the Cbs of 
hadrons reaches one. QCM predicts the position at fiB ~ 0.3 
GeV corresponding to the AA collisions at collisional energy 
about 12 GeV while statistical model predicts the position at 
fiB ~ 0.4 GeV corresponding to the collisional energy about 9 
GeV. The Beam Energy Scan programme of STAR collabora¬ 
tion can distinguish these two different predictions. 


VI. SUMMARY AND DISCUSSION 

Hadronization describes the process of the formation of 
hadrons out of quarks and/or gluons. In the meanwhile, the 
correlation properties between conservative charges of the 
system change also due to the transformation of basic degrees 
of freedom of the system. Phenomenological models mecha¬ 
nisms of the hadronization should reproduce these changes of 
charge correlation properties during hadronization process for 
the system. In this paper, we have studied the change of the 
baryon-strangeness correlation caused by hadronization in the 
quark combination mechanism. We found that this correlation 
is a good quantity of studying the hadronization because of 
its sensitivity on the dynamics of production competition be¬ 
tween baryons and mesons. We calculated the correlation co- 
efficientCBs = -3((B5)-(B){5))/{5^) of initial hadrons pro¬ 
duced from the deconfined free quark system with = 1. 
The calculated Cg^ of initial hadrons under Poisson fluctua- 
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tion is about 0.65 at zero baryon chemical potential, which 
is consistent with the prediction of statistical model for the 
hadron resonance gas. Beyond Poisson statistics, we calcu¬ 
lated of initial hadrons by IdS The 

resulting at zero baryon chemical potential is about 0.78, 
comparable to the calculations of Lattice QCD in the phase 
boundary at Tc - 163 MeV. This suggests that the quark com¬ 
bination, as a phenomenological microscopic mechanism, is 
able to describe the change of conservative charge correlations 
at hadronization, revealing certain basic dynamics of the real¬ 
istic hadronization process. We also predicted the correlation 
coefficient of hadrons at different baryon chemical potentials 
and compared them with the existed calculations of statisti¬ 
cal method. These predictions are expected to be tested by 
the Beam Energy Scan experiment of STAR Collaboration at 
RHIC and/or by the future Lattice QCD calculations at non¬ 
zero chemical potentials. 

Finally we would like to comment the measurement of Cbs 
directly by the definition Eq. ([TJ either in realistic experiments 
or in model simulations using an event generator, which al¬ 
ready has lots of discussions in literatures. To measure the 


Cbs in general we usually have to select a finite accept win¬ 
dow, e.g. mid-rapidity region |y| < 0.5. Selecting a small 
window size will increase the Poisson statistical fluctuations 
while selecting a wide window size will inevitably involve the 
global conservation effect of baryon number and in particu¬ 
lar of the strangeness. Amending such a finite window effect 
is usually complex. These considerations just prompt us to 
calculate the Cbs of the system in Sec. V by the response of 
baryon number of the system with respect to the change of the 
system strangeness. This method is less relevant to the accept 
window size and thus maybe more closely reflect the intrinsic 
baryon-strangeness correlation of the hadron system. 
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